HARNACK INEQUALITY FOR NONDIVERGENT PARABOLIC OPERATORS 
ON RIEMANNIAN MANIFOLDS 



Abstract. We consider second-order linear parabolic operators in non-divergence form 
that are intrinsically denned on Riemannian manifolds. In the elliptic case, Cabre proved 
a global Krylov-Safonov Harnack inequality under the assumption that the sectional cur- 
vature of the underlying manifold is nonnegative. Later, Kim improved Cabre's result by 
replacing the curvature condition by a certain condition on the distance function. Assuming 
essentially the same condition introduced by Kim, we establish Krylov-Safonov Harnack 
inequality for nonnegative solutions of the non-divergent parabolic equation. This, in par- 
ticular, gives a new proof for Li-Yau Harnack inequality for positive solutions to the heat 
equation in a manifold with nonnegative Ricci curvature. 
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1 . Introduction and main results 

In this paper, we study Harnack inequalities for solutions of second-order parabolic 
equations of non-divergence type on Riemannian manifolds. Let (M,g) be a smooth, com- 
plete Riemannian manifold of dimension n. For x e M and ( e R, let A x<t be a positive 
definite symmetric endomorphism of T X M, where T X M is the tangent space of M at x. We 
denote (X, Y) := g(X, Y) and \X\ 2 := (X,X) and assume that 

(1) A\X\ 2 < {A XJ X, X) < A\X\ 2 , V(jc, t) € M X R, Vie T X M 

for some positive constants A and A. We consider a second-order, linear, uniformly para- 
bolic operator _§f defined by 

(2) Jz?m = Lu-u, :- trace(A v o D 2 u) - u, in MxE, 

where o denotes composition of endomorphisms and D 2 u denotes the Hessian of the func- 
tion // defined by 

D 2 u ■ X = V x Vm, 

where Vu(x) e T X M is the gradient of u at x. Notice that in the special case when A x4 = Id, 
the equation ££ u = simply becomes the usual heat equation u, - Am = 0. 

In the elliptic setting, Cabre proved in a remarkable paper llCal that if the underly- 
ing manifold M has nonnegative sectional curvature, then Krylov-Safonov type (elliptic) 
Harnack inequality holds for solutions of uniformly elliptic equations in non-divergence 
form. Later, Kim (K| improved Cabre's result by removing the sectional curvature as- 
sumption and imposing a certain condition on distance function which, in the parabolic 
setting, should read as follows: For all p e M, we have 

(3) Ad p (x) < H f^- for x t Cut(p) U {p}, 

d p (x) 

Or 

(4) LdJx) < — — for x $ Cut(p) U \p), t e R, 

d p (x) 

where d p (x) = d(p, x) is the geodesic distance between p and x, Caiip) denotes the cut 
locus of p, and cil is some positive constant that is fixed by the operator L. We shall 
prove that if the above conditions ([3]) and Q hold, then we have Krylov-Safonov Harnack 
inequality for the parabolic operator ££; i.e., if u is a (smooth) nonnegative solution of 
Jzf w = / in a cylinder Kir := Bydxo) x (fo - 4R 2 , to), where x e M and ?o 6 R, then we 
have 

(5) sup M <c{inf M + W / f 

k s \n \voi(/sr 2 fi) Jk 2R / J 

where K~ :- B R (xo) x (to - 3R 2 , t - 2R 2 ), := B R (x ) x (t - R 2 , to), Vol denotes the 
volume, and C is a uniform constant depending only on n, A, A and a^. It is well known 
that the condition |3} holds if the manifold M has nonnegative Ricci curvature. Also, as it 
is proved in [K|, the condition Q is satisfied, for example, if for all x € M and any unit 
vector e € T X M, we have Mr \R(e)\ > 0. Here, R(e) is the Ricci transformation of T X M 
into itself given by R{e)X := R(X, <?)<?, where R{X, Y)Z is the Riemannian curvature tensor, 
and 

M~[R(e)] = M~[R(e), A, A] := A Y k { + A Y ^, 

where a:, are eigenvalues of the (symmetric) endomorphism R(e). In the case when Jz? is 
the heat operator and M has nonnegative Ricci curvature, then the condition Ai~[R(e)] > 
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is satisfied and thus the Harnack inequality |5]l holds; i.e., if M has nonnegative Ricci 
curvature, then we have 

sup M < cJinf u + R 2 ( f \u, - Au\ n+l \ \ , 

k« \VaKK21d J J ) J 

where C„ is a constant that depends only on the dimension n. This, in particular implies 
the Harnack inequality of Li and Yau HLY1 Also, in the case when M has nonnegative 
sectional curvature, then the condition Ai~[R(e)] > is trivially satisfied and we have the 
inequality |5]) with a constant C depending only on n, A, A, which especially reproduces 
the Harnack inequality by Krylov and Safonov [KS | in the Euclidean space. 

One crucial ingredient in proving the Euclidean Krylov-Safonov Harnack inequality is 
the Krylov-Tso estimate, which is the parabolic counterpart of the Aleksandrov-Bakelman- 
Pucci (ABP) estimate. The Krylov-Tso estimate as well as the classical ABP estimate is 
proved using affine functions, which have no intrinsic interpretation in general Riemannian 
manifolds. In the elliptic case, Cabre ingeniously overcame this difficulty by replacing the 
affine functions by quadratic functions; quadratic functions have geometric meaning as the 
square of distance functions. Following Cabre's approach, we introduce an intrinsically 
geometric version of Krylov-Tso normal map, namely, 

3>(x, f) := (exp t V x u(x, t),——d(x, exp t Vm(x, f)) 2 - u(x, t) 

The map <t> is called the parabolic normal map related to u(x, t). A few remarks are in 
order regarding the normal map. In the classical ABP (and Krylov-Tso) estimate, an affine 
function concerning with the (elliptic) normal map x h-> Vu(x) plays a role to bound the 
maximum of u by estimating the measure of the image of the normal map. Since an affine 
function cannot be generalized naturally to an intrinsic object in Riemannian manifolds, 
Cabre used paraboloids instead in [Ca|. The map 

p i-> min{M(x) - p ■ x] for a domain Q 

is considered (up to a sign) as the Legendre transform of u. Krylov |Rr] discovered the 
parabolic version of the Aleksandrov-Bakelman maximum principle and Tso [T| later sim- 
plified his proof by using the map 

(x, t) i — > (V x u(x, t), V x u(x, t ) ■ x — u(x, /)). 

We end the introduction by stating our main theorems. The rest of the paper shall be 
devoted to their proof. Below and hereafter, we denote 

f Q f ' = Vol(0 J/ 

and 

K r (x , t ) := B r (xo) x (t - r 2 , f ], (x , t Q ) e M x R. 

Theorem 1.1 (Harnack inequality). Suppose conditions (|3j, Q hold. Let u be a nonnega- 
tive smooth function in K2r{xo,4R 2 ), where Xq 6 M and R > 0. Then, we have 

(6) sup u<c\ inf u+R 2 (-f \^u\" +l ) I , 

K R (x ,2R 2 ) \K R (xoAR 2 ) \Jk 2 r(.^AR 2 ) I j 

where C is a uniform constant depending only on n, A, A and a^. 
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Theorem 1.2 (Weak Harnack inequality). Suppose the conditions (|3), Q hold. Let u be a 
nonnegative smooth function satisfying Jzfii < / in K2r(xq,4R 2 ), where xq e M and R > 0. 
Then, we have 

If A <c{ inf ^u+R 2 lf |/T +1 )" +I j; / + :=max(/,0), 

\Jk r (x ,2R 2 ) I {K R (x ,4R 2 ) \Jk 2 r(xoAR 2 ) 1 J 

where p e (0, 1) and C are uniform constants depending only on n, A, A and a^. 

2. Preliminaries 

Let (M, g) be a smooth, complete Riemannian manifold of dimension n, where g is the 
Riemannian metric and Vol := VoL, is the reference measure on M. We denote (X, Y) :- 
g(X, Y) and \X\ 2 := (X,X) for X, Y e T X M, where T X M is the tangent space at x e M. Let 
d(-, •) be the distance function on M. For a given point y e M, d y (x) denotes the distance 
function fromy, i.e., d y (x) := d(x,y). 

We recall the exponential map exp : TM — ■> M. If y Xr x : M. — > M is the geodesic starting 
from x e M with velocity X e r,M, then the exponential map is defined by 

exp,(X) := y XiX {\). 

We note that the geodesic y X) x is defined for all time since M is complete. Given two 
points x,y € M, there exists a unique minimizing geodesic exp K {tX) joining x to y with 
y = exp x (X) and we will write X = exp^ 1 ^). 
For X e T X M with |Z| = 1, we define 

t c (X) :- sup [t > : exp x (sX) is minimizing between x and exp x (tX)} . 

If f c (X) < +oo, exp J .(f c (X)X) is a cut point of x. The cut locus of x is defined as the set of 
all cut points of x, that is, 

Cut(x) := {exp x (t c (X)X) : X e T X M with \X\ = 1, f c (X) < +oo) . 

Define 

E x := {fX e T X M : < t < t c (X), X e r r M with |X| = 1} c T X M. 
One can show that Cut(x) = exp x (dE x ), M = exp^ZS j) U Cut(x) and exp v : E x — > exp^E^) 
is a diffeomorphism. We note that Cut(;t:) is closed and has measure zero. For any x t 
Cut(y) with x + y, then d y is smooth at x and the Gauss lemma implies that 

„ , , , exp-'(y) 

WdJX) = ; 

} |exp-i(y)| 

and 

V(4/2)W = -exp; 1 (y). 

Let the Riemannain curvature tensor be defined by 

R(X, Y)Z = V z V y Z - V y V x Z - V KF] Z, 

where V stands for the Levi-Civita connection. For a unit vector e e T X M, R(e) will denote 
the Ricci transform of T X M into itself given by R(e)X := R(X, e)e. 

For u € C°°(M), the Hessian operator D 2 u(x) : T X M -> T^M is defined by 

D 2 m(jc) • X = V x Vm(x). 

Let M and jV be Riemannian manifolds of dimension n and cp : M — » be smooth. The 
Jacobian of <p is the absolute value of determinant of the differential d<p, i.e., 

Jac (p(x) := | det d<f>(x)\ for xeM. 



HARNACK INEQUALITY FOR NONDIVERGENT PARABOLIC OPERATORS 



5 



We quote the following lemma from Lemma 3.2 in MCal . in which the Jacobian of the map 
x h-> exp x (Vv(x)) is computed explicitly. 

Lemma 2.1 (Cabre). Let v be a smooth function in an open set Q of M. Define the map 
4> : D. M by 

(f>(x) := exp c Vv(x). 
Let x € Q. and suppose that Vv(x) G E x . Set y :— <p(x). Then we have 

Jac <p(x) = Jac exp c (Vv(jt:)) ■ | det D 2 (v + d 2 /2j (x) |, 

where Jac exp v (Vv(x)) denotes the Jacobian q/exp v , a map from T X M to M, at the point 
VvU) e T X M. 

Under the condition ([3]), we have the estimate for Jacobian of the exponential map and 
Bishop's volume comparison theorem as follows. We state the known results as a lemma. 
The proof can be found in [K p. 286] (see also flQ). 

Lemma 2.2. Suppose that M satisfies ([3]). 
( ij For any x e M and X € E x , 

Jacexp x (X) = |detafexp x (X)| < 1. 

(ii) (Bishop) For any x € M, Vol(Z?«(x))//?" is nonincreasing with respect to R, where 
Br(x) is a geodesic ball of radius R centered at x. Namely, 

< — if < r < R. 

Vol(B r (x)) ~ r n J 

In particular, M satisfies the volume doubling property; i.e., Vol(Z?2«(*)) ^ 2" WoI(Br(x)). 

The following is the area formula, which follows easily from the area formula in Eu- 
clidean space and a partition of unity. 

Lemma 2.3 (Area formula). For any smooth function <p : M x R — > M x R and any 
measurable set E C M X R, we have 

\ Jac <p{x, t)dV(x, 0=| 'H\E n <f\y, s)]dV(y, s), 

JE JmxR 

where f H° is the counting measure. 

Notation. Let us summarize the notations and definitions that will be used. 

• Let r > 0,p > 0, z e M and t a e R. We denote 

K r , p (z , Q := B r (z ) x (t - p, t ], 

where B r (z ) is a geodesic ball of radius r centered at z - 

• We denote K r (z , t ) := K r r i{z , t ). 

• We say that a constant C is uniform ifC depends only on n, A, A and a L . 

• We denote f Q f:= ^ f Q f. 

• We denote \Q\ := Vol(Q). 

• We denote the trace by tr. 
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3. Key lemma 

In this section, we obtain Aleksandrov-Bakelman-Pucci-Krylov-Tso type estimate (Lemma 



3.2 1 for parabolic Harnack inequalities. We begin with direct computation of the Jacobian 



of the parabolic normal map below, which is a parabolic analogue of Lemma 2.1 



Lemma 3.1. Let v be a smooth function in an open set K of M x R. Define the map 
4>:K^Mby 

<f>(x, t) : = exp x V x v(x, f) 

and the map O : K -> M xRby 

<£>(x, t) := i^>(x, t),-^d (x, (f>(x, t)) 2 - v(x, t) 

Let (x, t) e K and assume that V x v(x, t) e E x . Set y :— cf>(x, t). Then 

Jac 0>(x, t) = Jac exp v (V A v(x, f)) • |(-v f ) det {pi (v + dy/l)) |, 
where Jac exp x (V x v{x, t)) denotes the Jacobian ofexp x at the point V ' x v{x, t) € T X M. 

Proof. We may assume that V ' x v{x, f) + 0, which is equivalent to x + y. Let (£, cr) e 
T X M x R\{(0, 0)} and let y = {y\ , y 2 ) be the geodesic with y(0) = (x, t) and y'(0) = (£, cr). 
We note that j\(f) - exp t t£ and y 2 (r) - t + err. Set 

F(i,r):=exp 7i(T) [iV x v( r (T))]. 

Consider the family of geodesies (in the parameter s) 

n(s,r) := (y(s,t), 72 (t) - sl^d(ji{T),cf>(y(T))) 2 + v(y(r)) + y 2 (r) 

that joins 11(0, r) = y(r) to 11(1, t) = <t>(y(T)). Then we define 

J{s):=—\ Tl{s,T), 

which is a Jacobi field along 

X(s) :- {&xp x (sV x v{x, t)), t - s <Xd{x, (f>(x, f)) 2 + v(x, t) + t 
Simple computation says that 

7(0) = (£cr) and 7(1) = — Myir)) = d®(x,t) ■ (£cr). 

OTIt=Q 

We also have 

DJ(0) = (D 2 v(x, f) £ + crV x v t (x, t) , -oy t (x, - cr 

- (V, (d*/2) (y), dexp,(V,v(jc, 0) • (d* (v + d*/2) (jc, f) • f + crV,v,(x, f))» . 
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In fact, we have 

dill dlli 
DJ(0) = D s — \ = D T — \ 

8t Ij=o,t=o ds U=o,t=o 

= D t | t=q (v,v(y(T)), - irf(yi(T), 0(y(r))) 2 - v(y(r)) - y 2 (T)) 
= (z^vCje,0-? + ff-V,v ( (x,0, 

- (V. v (4/2)(x), f) - h x (d 2 x /2)(<f>(x, 0), ^(T(T))| T= J - <V x v(x, f), f> - crv, - <r) 
D 2 t v • f + crV. c v„ - /v,(^/2)(0(x, f)), §^<K7(t))\ t 



to obtain 



since V x ((B/2)(x) = - exp x 1 (y) = — V x v(x, ?)■ Then we use Lemma 2. 1 

— 0(rW)| T _ o = </exp x (V,v(x, f)) • (D* (v + 4/2) (x, f) ■ f + <rV,v,(x, f)) ■ 
On the other hand, consider the Jacobi field along X(s) satisfying 
7 fj(r (0) = (£cr) and 7^(1) = (0,0). 

Then we can check that 

/f^(j) = ^^| t=0 and D s 7 fi£r (0) = (4/ 2 ) (*) ' 6 ~°) • 

where 

T(i, t) := ^exp yi(T) s exp" 1 ^ f), Ji(t) - s l~d(x, <f>{x, t)f + v(x, t) + y 2 (r) jj . 

(We refer [Ca, Lemma 3.2] for the proof.) 

Define := 7 — J^g-. The Jacobi field 7f j(r along X(s) satisfying 

= (0- 0) and 0,7^(0) = D s 7(0) - D s J^(0) 

is written by 

dexp (xl) (sX'(0))-(sDj^(0)). 

Therefore, we have 

7(1) = 7^(1) = dap M h x v(x, t),- l -d(x,yf - v(x, t) - t J • (0,7(0) - 0,7^(0)) ., 
which means 

d<£>(x, t) ■ (£ &) = d exp (j . t) (v x v(x, f), - ^d(x, y) 2 - v(x, t) - t 
(D 2 x (v + d 2 /2)(x,t)^ + (rV x v t (x,t), 

-(tv, - (V, (d 2 J2) (y), dezp x (V x v(x, t)) ■ (d 2 x (v + d 2 /2) (x, t) ■ arV x v t (x, f)))) 

= (dexp x (V x v(x, 0) • (D 2 X (v + d 2 /2) (x, t) ■ £ + aV x v t (x, tj) , 

-oy t - (V, (d 2 J2) (y), dsxp x (V x v(x, t)) ■ (d 2 x (v + d 2 /2) (x, t ) • f + <rV,v t (*, f)))) . 



s 
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To calculate the Jacobian of <t>, we introduce an orthonormal basis \e\, ■ ■ ■ , e„) of T X M and 
an orthonormal basis {ei, ■ • ■ , e„} of T V M = T ex? ^ v ^ xt )M. By setting for i,j= 1, • • ■ ,n, 

A tJ := («,, dexp x (V x v(x, t)) ■ (d 2 x (y + d 2 /2) (x, t)ej)) , 
b/ : = (e,-, dexp x (V x v(x, ?)) ■ Vv t (x, t)) , 
Ci :=(e h V x (d 2 x /2)(y)), 
the Jacobian matrix of <E> at (x, f) is 

Ay bi 

-CkA k j -v,-b k c k 
Lastly, we use the row operations to deduce that 

An bi 



Jac ®(x, 
This completes the proof. 



det 



-v, 



-v f )det(Ay)|. 



The following lemma will play a key role to estimate sublevel sets of u in Lemma |4.3| 



and then to prove a decay estimate of the distribution function of u in Lemma 6.1 This 
ABP-type lemma corresponds to ICal Lemma 4.1]. 

Lemma 3.2. Suppose that M satisfies the condition Let z & M, R > 0, and < r/ < 1. 
Let u be a smooth function in K aiRa2R t(z , 0) C M X R satisfying 

(7) m>0 in K aiR ^ 2R 2(z o ,0)\Kfj lR ^ 2R 2(zo,0) and inf h < 1, 

where a\ :— —, a,i : = 4 + rj 1 + \, B\ :— -, and B2 '■= 4 + rj 2 . Then we have 



(8) 



\B R (zo)\ R 1 < C(n,n,A) \ \{r 2 %u + a L + A + l) + 

J|«<M„)n^ iSfcs2 (z„,0) I 



re+1 



where the constant M n > depends only on r\ > one/ C(/7, n, A) > depends only on rj, n 
and A. 

Proof. For any y e B R (z ), we define 

1 2 1 2 6 

Wy(.x, f) := -fl u(x, t) + -d-(x) - C v t, C n := -r. 

2 2 - if- 



(zo,0) 



2-R 





1 


4i? 2 






/,-'/,-' 





a 2 7? 2 



Figure 1. ai := ^,a 2 := 4 + rf + \,Bi := 9 -,B 2 := 4 + i] 2 
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From the assumption Q, it is easy to check that 

inf Wv < |5 + ^t\R 2 --A„R 2 , 

K2r(Zo,0) \ Tf j ' 

and 

Wy > ^6 + R 2 = (A n + l)R 2 on K aiKa ^{z B , 0)\K m/h Az o , 0). 

From the above observation, for any (y, h) € B R (z a ) x ^A^ 2 ,(A 7 + l)^ 2 ), we can find a 
time t e (-B 2 R 2 ,Q) such that 

h = inf wy(z, t) = Wy (x, t] , 

where the infimum is achieved at an interior point x of Bp lR (z )- By the same argument as 
in BCal pp. 637-638], we have the following relation: 



y = exp x V^fl^J (x~i). 



Now, we consider the map 4> : K aiR ^ R i(z , 0)^MxR (with v(jc, f) = ^R 2 u(x, t) - C n t 
in Lemma |3~T) defined as 



OO, f) : = 
Define a set 



1_, \. . 1 ./ _ /!_, \. ,\ 2 1 



exp t V A \-R u Ux, t),--d\x,e\p x V v \-R u Ux, t)\ - -R u(x, t) + C n t 



E:=Ux,t)€ K PlR ^ R i{z o ,0) : 3y e B R (z„) s.t. w y (x, t) = inf w y < (A, + \)R~ 

The set £ is a subset of the contact set in Kp xR ^ R i(z , 0) that contains a point (x, f), where 
a concave paraboloid -^d 2 (x) + C n t + C (for some C) touches \R 2 u from below. Thus 
we have proved that for any (y, s) e B R (z ) x (-(A ;; + l)R 2 , -A n R 2 ), there is at least one 
(x, f) e E such that (y, s) = 0(x, f), namely, 

Br(Zo) x (-(A, + l)fl\ -A n R 2 ) c <D(£). 

So Area formula gives 

(9) \B R (z )\ R 2 < f H \E n _1 Cy, s)| dV(y, s) = f Jac 0>(x, 0dV(jc, f). 

We notice that for (x, t) e Eandy e B R (z ), w y (x, t) = \R 2 u(x,i)+\d 2 {x)-C ri t < (A,+l)/? 2 
and hence u(x, t) < 2{A n + 1) =: M v for (x, t) e E. 
Lastly, we claim that for (x, t) e E, 



(10) Jac 0(x, r) < - - \^-R 2 J?u(x, t) + a L + A + C, 

Fix (x, t) e E and y e B R (z ) to satisfy 

w y (x,0 = inf Wy. 

We recall thaty = exp x V, {\R 2 u) (x, t) (see [Ca, pp. 637-638]). 
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If x is not a cut point of y, then Lemma 3.1 (with v(x, t) = ^R 2 u(x, t) - C^t) and Lemma 



2.2 (i) imply that 

Jac 0(x, f) < 



1 



R-it, + C„ ) det|z) 2 | -A'-// + -<l; j ) (.v. / 1 



1 



1 



Since the minimum of w y in Bp^iz,,) x {-B2R 2 , t] is achieved at (x, t), we have 

< D 2 t w y (x, t) = D\ {^R 2 u + and - 5 ? w y(*> f ) = ^ 2 "' _ C >7> 

where D 2 w y (x, f) > means that the Hessian of w y at (x, f) is positive semidefinite. There- 
fore, by using the geometric and arithmetic means inequality, we get 

Jac (D(x, < {^-^R 2 u, + C,j det [d 2 x + -dyjj (x, f) 

^ 2 M , + C,JdetA v detjD 2 ^/? 2 M + ^</ 2 



1 / 1 



A" \ 2 
1 



(n + l) n+1 A" 
1 

(n + 1)"+U" [2 



1 , 1 „ 



— tr A,, o -R>u + - -Jtti, + C, 



1 



R 2 ^fu(x, t) + L 



2 y 



1 



2 



/? 2 Jz^m + a L + A + C, 



where we used 



(n + 1)" + U" [2 

= 1 — r - I [\r 2 ^u + a L + A + C„) 

(n + l)" +l A n \\2 n ) 



L [d 2 /2] = d y Ld y + (A XJ Vd y , Vd y ) <a L + A|VJ V | 2 . 



When x is a cut point of y, we make use of upper barrier technique due to Calabi |Cal|. 
Sincey = exp v V x (jR 2 uj (x, t), x is not acut point ofy^ := (f>o-(x,t) '■= sxp x 'V x ^jR 2 uj(x,t) 
for < <x < 1 . Now we consider 



®Az,t) := \4> <t {z,t),--R l u(z,t)- -d(z,Mz,T)Y + C n T 

instead of since Jac i>(x, t) = lim Jac <t>, r (x, t). As before, we have 

o-Ti 



Jac fl&o-Cx, t) < 



(-^R 2 u, + C, ; )det(z) 2 (^ 2 M + l -d], 



(x,t) 



We note that 



lim inf 

o-Tl 



(-^R 2 u, + C,)det|z) 2 |^ 2 M + l -d\ 
= liminf \[-d,w y Jx, f)) det (D 2 x w yir ) (x, f)| 



(x,t) 
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for w y<r (z, t) : = \R 2 u{z, t) + \d 2 ytr {z) - C n r. According to the triangle inequality, we have 

w y (z, t) < l -R 2 u{z, (dyjz) + d(y a ,y)f - C n r 

1 2 

= w >v (z,t) + d(yo-,y)dy a (z) + -diy^y) , 

where the equality holds at (z, r) = (x, t). Since w y has the minimum at (x, f) in B^r{z ) X 
(-P2R 2 , t], the minimum of iv fr (z, T)+d(y cr , y)d y<T (z) (in B / 3 1 r(z )x(-;S2^ 2 , ?]) is also achieved 
at (x, t), that implies that 

D 2 {w yiT + (iOo-, y)c/ V(r ) (x, > 0, dtWyJx, t) < 0. 

To bound D 2 y cr (x) uniformly in cr e [1/2, 1), we recall the Hessian comparison theorem 
(see llSl USYIQ : Let -k 2 (k > 0) be a lower bound of sectional curvature along the minimal 
geodesic joining x and y. Then for < cr < 1, 

D 2 d y Jx) < kcoHh(kd y<r (x))Id 

and hence we find a constant N > independent of cr such that 

D 2 d y<r (x) < A/Id for - < cr < 1. 

Following the above argument, for j < cr < 1 , we obtain 

< lim inf (-d t w y (x, t)) det (D 2 w y + dOv, y)D 2 d y ) (x, f) 

< lim inf (-d,w >v (x, f)) det (d 2 w, v + dtja-, y)Md) (x, f) 

1(1 1 " +1 

< lim inf : — \ -R 2 J£u + a L + A + C„ + d(yg-,y)nAN } 

<ni (n + \2 ' J 

„ , 1 

— — \\-R 2 ^u + a L + A + C„ | 



" ( 

Then we deduce that 



Jac <£(x, t) < — ; — U-R 2 £fu(x, t) + a L + A + C„ | 







lim inf det (d 2 w y ) (x, f) = lim inf det (d 2 w v + d(ya-, y)Nld) (x. 

= lim mf det (£> 2 vv >v + d(y a , y)N\&) (x, t). 

We conclude that (TO) is true for (x, f) e £\ Therefore the estimate ([8) follows from (|9]l 
since £ c \u < M^jn ^ lS , A «2(z o ,0). □ 



4. Barrier functions 



We modify the barrier function of [ W | to construct a barrier function in the Riemannian 
case. First, we fix some constants that will be used frequently (see Figure [TJ; for a given 
< 77 < 1, 

11 , 77 4 9 , 

a\ := — , a 2 := 4 + 77 H , j8i := - and y8 2 := 4 + 77 . 

77 4 77 
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Lemma 4.1. Suppose that M satisfies the condition Q. Let z n e M, R > and < 

rj < 1. There exists a continuous function v n (x, t) in K a]R a2R 2(z ,(3 2 R 2 ), which is smooth in 
(M\ Cut(z )) n K aiRta2R i(z ,/3 2 R 2 ) such that 

(i) v„(x, t) > in K a ^ aiR 2{ Zo ,p 2 R 2 ) \ Kp 

iR,/3 2 R 2 

(ii) v^x, t) < in K 2R (z ,f3 2 R 2 ), 

(Hi) R 2 ^v„ + a L + A + 1 < a.e. in K 0lR ^ lR 2(z o ,/3 2 R 2 )\Ki R (z o , ^R 2 ), 

(iv) R 2 ^ Vll < C„ a.e. in K PlR , hR 2(z ,fi 2 R 2 ), 

(v) v n (x, t) > -C ;; in K aiR a2R i(zo,p 2 R 2 ). 

Here, the constant C n > depends only on T], n, A, A, a L ( independent ofR and z )■ 

Proof. Fix < r] < 1 . Consider 

1 



h(s,t) :- -Ae 



1 



fx) 



exp 



K) 



for t > 0, 



(47Tf)" /2 

as in Lemma 3.22 of [W| and define 

ifr(s, t) := h(s, t) + (a L + A+l)t in [0,#] x [0,J3 2 ]\[0, £ ] X [0, £], 

where the positive constants A,m,l,a ( depending only on T],n,A,A,a L ) will be chosen 
later. In particular, / will be an odd number in N. We extend iff smoothly in [0, ffj]x[- -4- , jS 2 \ 
to satisfy 

*A > on [0,a*| x [-f.&MO,^] x [0,&], 
<A>-C„ on [0, oft x 



and 



sup {2a L + A (2|5 S <A| + 4j|3„fl) + |<9^|} (j, t) < C, 

[0,^]x[0j3 2 ] 



for some > 0. We also assume that iff(s, t) is nondecreasing with respect to s in [0, a 2 ] x 
[—j-,f> 2 ]. We define 



l d 2 (x) t 



v n (x, t) - v(x, t) :— 



for (x, f) € K aiR a „g2 



R 2 'R 2 t 

where d Za is the distance function to z a . Properties (i) and (v) are trivial. 

We denote d u (x) and h ^ , j by d(x) and <f>(x, t) for simplicity and we notice that 

for (x, t) e K PiR „ 2R 2(zo,p 2 R 2 )\K u (zo, ^R 2 ), 



v(x, t) = h 



d 2 (x) t 
~R 2 ~ , ~R 2 



+ (a L + A + 1)— = </>(x, t) + (a L + A+ 1)^ 



and (p{x, f) is negative in Kp [R p 2R 2 
Now, we claim that 



(11) 



< a.e. in K m ^ R 2(z ,/3 2 R z )\K lR (z , \R l ). 

Once ( fTT| is proved, then property (iii) follows from the simple calculation that R 2 J/? [(aL + A • 
-(a L + A + 1) in Kp lR p 2R 2(zo,/3 2 R 2 ). Now we use the identity 

Jt?[ip(u(x), t)] = d u ip(u, t)££u + d uu ip(u, t){A XJ Vu, Vm> - d,tp(u, t) 
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to obtain 

2d ^ Id 2 t\ 

2 4d 2 1 Id 2 t \ 1 id 2 t 

-,d s h + —rd ss h \ [ T , T I (AwW, Vd) - —d t h 
R 2 R 4 )\R 2 R 2 j ' R 2 \R 2 R 2 J 

Since d ■ Ld < a L and A < (A Xit Vd, Vd) < A in M\ Cut(z ), we have that 

Q&R 2 _ d 2 ) 2 



(-4 



2a 1 
)— h 



= (ft/? 2 - d 1 ) |2/ + (ft/? 2 d££d 

- \l(l- 1)4</ 2 + 2l(/3 2 R 2 -d 2 )^- + (ft 2 /? 2 -d 2 ) 2 ^- \<A u Yd.Vd) 

+ (ft 2 /? 2 - d 2 ) ^21 + (J3\R 2 - d 2 ) y j (A x>t Vd, Vd) 
, ad 



+ <J3 2 R 2 -d 2 ) 2 ^L-<fi 2 R 2 -d 2 ) 2 (^ + ^) 



i 2d' cxd 

< 2l(J3 2 R 2 - d 2 )(a L + A) + OS 2 /? 2 - d 2 ) 2 I —(a L + A) + — 

t i i to/ 4-a 2 d 2 n 

- 1(1 - l)4d 2 A - (ft 2 /? 2 - d 2 ) 2 — - — A + - 

\ t l 2t 

4cKt/^ hi 

- 2l(J3 2 R 2 - d 2 )—A - —(J3 2 R 2 - d 2 ) 2 a.e. in /^^(z^ft/? 2 ). 



By choosing 



(12) 



we deduce 



1 2£i 

a := — , -^-(a L + A) + 1 < / := 2/ + 1 (for some I' e N), 
4/1 77 2 /l 

„ J 8a, 2/( flL + A) | 
w := 2 • max<{ —t(cil + A), : — V, 



~ -J^0<O a.e. in^ lfisS2R 2(z ,yS 2 /? 2 )\Wz ,^ 2 ). 



(-0) 

' N /( A'./i; 



Indeed, we divide the domain Kp iR R^ R i(z , P2R 2 )\K>i R (z n , \R 2 ) into three regions such that 



„2 

K PxRfilR 2(xo,P2R 2 )\K iR (z , jR 2 ) =: Aj U A 2 u A 3 , 

where Ai := {0 < £ < f < f < ft}, A 2 := < £ < ft, § < f < ft } and 
A 3 :={^ < <ft,0< \ < §}. We can check that 

(J3 2 R 2 -d 2 ) 2 , „ 2 . 

^— ^0 < a.e. in /^ lRAs2 (z ,ft/? 2 )\/i: !R (z , ^/? 2 ) 

by choosing a and I large in A 1? m large in A 2 and A3 as in ( p"2| ). Therefore, we have proved 
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From the assumption on ifr, we have that for a.e. (x, t) e Ka lR a lR i{z , B2R ), 



R 2 ^v(x, t) = 2d s if/ 



d 2 t 



. Ad 2 
d ■ Ld + \ 2<9 s i/r + -jp-d ss ijf 



R 2 'R 2 



(A x , t Vd, Vd) 



< sup {2a L d s i// + A (2d s ifr + 4s\d ss iff\) + \d t iff\] (s, t) < C n . 

This proves property (iv). 

In order to show (ii), we take A > large enough so that for (x, t) e KwizofhR 2 ), 



v(x, t) < -Ae~ hm 



1 - 



1 



-4a I if 



(4tt/? 2 ) 



n/2" 



+ (a L + A+l) y 6 2 <0. 



This finishes the proof of the lemma. 



Now we apply Lemma 3.2 



to u + v n with v l; constructed in Lemma 4. 1 and translated 
in time. Since the barrier function v n {x, t) = ip n y*§r-, -gi j is not smooth on Cut(z a ), we 

need to approximate v, ; by a sequence of smooth functions as Cabre's approach at HCal . 
We recall that the cut locus of z is closed and has measure zero. It is not hard to verify the 
following lemma and we just refer to HCal Lemmas 5.3, 5.4. 

Lemma 4.2. Let z 6 M, R > and let tfi : R + x [0, T] — > R be a smooth function such 
that i/f(s, t) is nondecreasing with respect to s for any t € [0, T]. Let v(x, t) := tp (<? (x), f). 
Then there exist a smooth function < £(x) < 1 on M satisfying 

g = 1 in Ba R (z ) and supp£ c Bm R (z ) 

'i 

and a sequence {wk\™ =l of smooth functions in M X [0, T] such that 

Wk — > f v uniformly in M X [0, T], 

d t wu — * £d t v uniformly in M X [0, T], 

D 2 x w k < CId in M x [0, T], 

D 2 x w k -> D\v a.e. in B PiR ( Zo ) x [0, T], 

where the constant C > is independent ofk. 

Lemma 4.3. Suppose that M satisfies the conditions ([3]),Q. Let z G M,R > 0, and 

< t] < 1. Let u be a smooth function such that ££u < f in K a ^ a ^ R 2(z ,4R 2 ) such that 

u>0 in K aiRaiR 2 izoAR 1 ) 

and 

inf u < 1. 

Then, there exist uniform constants M„ > 1, < fl v < 1, and < s n < 1 such that 

\{u < M„} n K, lR (z„,0)\ 



(13) 
provided 

(14) 



k,v,R,a 2 ^fe,4^ 2 ) 



"If 



1/ 



+ IH+1 



< S„ 



^ lS ,„ 2 «2fe,4« 2 ) 

where M q > 0, < e n , ji n < 1 depend only on tj, n, A, A and a L . 
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Proof. Let v n be the barrier function in Lemma 4. 1 after translation in time (by -77 R ) 
and let {w^j^lj be a sequence of smooth functions approximating v n as in Lemma 



4.2 



We notice that u + v„ > in a „ R z(Zo, 4R z )\K BlR B , R i{z , 4R Z ) and inf (u + v„) < 1. 

Thanks to the uniform convergence of wu to £v n , we consider a sequence {e/tl^j converging 
to such that sup Wk < £* and 

K 2 r(z„AR 2 ) 

w k > -s k in K a ^ a2R 2{z AR 2 )\K l3l 

R,/3 2 R 2 (zoAR 2 ), 

and define 

_ U + Wk + Sk 

1 + 2s* 

Then satisfies the hypotheses of Lemma 3.2 (after translation in time by 4R 2 ). Now we 
replace u by Wk in ([8]l and then the uniform convergence implies that for a given < 5 < 1, 
we have 



\B R {z )\R 2 < C(r,. 



n, A) \ [(R 2 ^fw k + a L + A + 1 V 



if k is sufficiently large. Since D 2 Wk < Cld and \d,Wk\ < C uniformly in k on Kp iR ^ R 2(z AR 2 ), 
we use the dominated convergence theorem to let k go to +00. Letting 5 go to 0, we obtain 

\B R (z )\ ■ R 2 < C(J/, «, /I) f i(R 2 ^[u + Vn ]+a L + A+ 

J{«+v,<M,}nir ftw2 fc,4fi2) I > 

>n+l 



C(j],n,A) \ l(R 2 ^[u + v n ]+a L + A+l) + 

j£,U£ 2 ^ 



where := {u + v n < M,,} n (^ lS>AS 2(z o ,4tf 2 )\^«(z o ,0)) and £ 2 := {m + v, < M ;; ) n 
K nR (z , 0). From properties (iii) and (iv) of \ n in Lemma [4~T] and Bishop's volume compar- 
ison theorem in Lemma 12. 21 we deduce that 



\K aM A, AR 2 )^ ± C V ||(^«rL (w(w(1(1)) + c ? lk*IU( wfc ^) 

where C ); > depends only on n, A and 77 > 0. We note that \u < M n - v n ] c \u < M ); + C n ) 
from (v) in Lemma 4. 1 Therefore, by taking 

1 JL 1 

S " = 2C~' M ^ = M " + C " and ^' + ' = 2C~' 
|{M<M;}n^«( Zo ,0)| 

we conclude that — ; ; — > u„ > 0. □ 

\K mR , a2R 2(z AR 2 )\ 



Using iteration of Lemma 4.3 we have the following corollaries. 
Corollary 4.4. Suppose that M satisfies the conditions ((3j, <|4]). Let z £ M and < 77 < 1. 

— /tV-I - ; —2 

For i e N, let R, :— y-J R and ti :— 2j}=i 4Rj. Let u be a nonnegative smooth function 
such that J£u < f in I ) k : _, K - -2(z„,f;) for some k e N. We assume that for h > 0, 
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inf u < h and 

Uii K 2lt XzJd 



■ If irr 1 



< s n hM k ~\ VI < ; < k. 



Then we have 
(15) 



\{u < hM k n ) n K nR ( Zo ,0)\ 

\K aiR , a2R 2(ZoAR 2 )\ 



■Mr,, 



where M n , £,,,//,, are the same uniform constants as in Lemma 4.3 

Proof. We may assume h = 1 since v := | satisfies Jzfv = ^-S?m < {■ We use the induction 
on k to show the lemma. When k = 1, it is immediate from Lemma [43] 

Now suppose that (JT3J is true for k - 1 . By assumption, we find a j e N such that 
1 < j„ ^ & and inf u — inf u < 1. Define v := ujMt"^ . Then v satisfies that 



< f/M; 



k-j, 
1 



inf v < 1 and 



\riM', 



T 2 <.Zo,tj ) 



< s„ 



Applying Lemma 



4.3 



to v in K - -i (z„,t; ), we deduce 



(v < M„ 


nK ,R,J^h-4Rl) 




(v < M, 


n ^fi fc _,Cz».*y.-i)| 

















which implies that inf u < inf u <M\ 
tion hypothesis for y' - 1(< — 1) to conclude 



: > 



Therefore, we use the indue - 



ulM k - h+x <M ] n ° 'jn^tzJ)) 
|^ff,fi,ff,« 2 fe,4^ 2 )| 



: yu, > 0, 



which implies ( |T5| . 



We remark that Lemma 4.3 and Corollary 4.4 hold for any M' n > M n . The following is 
a simple technical lemma that will be used in the proof of Proposition |4^6] 

Lemma 4.5. Let A,D>0 and s > 0. Let u be a nonnegative smooth function such that 
£fu < f in B R (z ) x (-AR 2 , 0] with 



R2 (f 



\f 



+ IH+1 



< s. 



„)x(-AR 2 ,0] 

Then, there exists a sequence Uk of nonnegative smooth functions in B R (z ) X (-AR 2 , DR 2 ] 
such that Uk converges to u locally uniformly in B R (z„) X (-AR 2 ,Q] and ££uk < gk in 
B R (z )x(-AR 2 ,DR 2 ] with 



R 2 lf 

\Jb, 



>B r (z )x(-AR 2 ,DR 2 ] 



\Sk\ 



< E. 
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Proof. First, we define for (x, t) e B R {z ) x (-00, DR 2 ], 

( for t € (-00, -AR 2 ], 

u(x,t):=< u(x,t) for te(-AR 2 ,Q], 

[ u(x,0) + St for te(0,DR 2 ], 

where S := sup {(Jz? u) + (x, 0) + \u t (x, 0)|}. Then u is Lipschitz continuous with respect to 

Br(Zo) 

time in B R {z ) x (-A7? 2 , DR 2 ] and satisfies 

C for r e (-00, -AR 2 ), 

Jfu(x,t) <f(x,t) := \ f(x,i) forf e(-AR 2 ,0), 

[ JS?k(jc, 0) + m,(jc, 0) - S < for f e (0, Dfi 2 ]. 

Let e k > converge to as £ — > +00, and let if be a nonnegative smooth function such that 
co(0 = for t g (0, 1) and £ ct>(f)<if = 1. We define coi(0 := — q?( — | and 

u k (x,t):- J u(x,s)<p k (t - s)ds, V(x, f) e fi R (z ) x (-00, D/? 2 ], 

where we notice that the above integral is calculated over (f - e^., f) c R. Then, a smooth 
function u k satisfies 

5£u k (x, t) = I ££u{x, s)tp k (t - s)ds < g k (x, ?), V(x, f) e B R (z a ) x (-00, DR 2 ], 



where g k (x, t) :- \ f (x, s)tp k (t - s)ds > 0. We also have 

R 2 



\Jb r (z,,)x(- 



AR 2 ,DR 2 ] 



\8t 



\B R ( Zo )\-(A + D)R 2 }-+ 
R 2 

\B R ( Zo )\-(A + D)R 2 }^ 

1 

A 



: 1 1/ I (B R (z„)x(-AR 2 -s k ,DR 2 -s k ]) 
fll/ llL"+'(B fi fc,)x(-AR 2 ,0]) 



U + D/ 



E < S, 



which finishes the proof. □ 

Proposition 4.6. Suppose that M satisfies the conditions ([3]),(|4]i. Let z s M, R > 0,0 < 

77 < j awif t e [3, 16]. Lef u be a nonnegative smooth function such that ££u < f in 
B^ R (z ) x (-3T? 2 , . Assume that inf m < 1 and 



JB«„fc,)x(-3i? 2 ,^l 
V IS v 1- J 



irr 



for a uniform constant < e'^ < 1. Let r > satisfy (^j R < r < f |) ^/or iVeN 
one/ /ef (zi, fi) fee a po/nf f/zaf t/(z ,Zi) < R and \t\\ < R 2 . Then there exists a uniform 
constant M' n > 1 (independent of r,N,Z\ and t\) such that 



u<M' N+2 }nK nr (zi,ti)\ 



\K„ 



2(zi,h +4r 2 ) 



: Hr, > 0, 
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where < fx n < 1 is the constant in Lemma 4.3 



Proof, (i) From Lemma 4.5 we approximate u by nonnegative smooth functions Uk, which 
are defined on Bm R {z ) x (-3T? 2 , (4 6 ^ 6 + 7? 2 ] . We find functions Uk and gk such that ut, 

converges locally uniformly to u in Bi» R (z ) x (-3R 2 , , and satisfies 



^u k < g k in B4| R (z ) x -3R 



64-R 2 



and 



fl 2 



\8t\ n+1 



(4 - ^ 2 )/7 6 
49 



< — e'< 2s' 
48 n " 



by using the volume comparison theorem and Lemma 4.5 For a small 6 > 0, we consider 
Wk '■= - ; and then for large k, Wk satisfies inf Wk < l,J^Wk < gk in Bis R (z )x 



1+6 



n+l 



according to the local uniform convergence of u k to u in Lemma 4.5 So if we show the 
proposition for Wk, the local uniform convergence will imply that the result holds for u by 
letting k — > +oo and 6 — > 0. Now we assume that u is a nonnegative smooth function in 
Bas r (z ) x (-3/? 2 , (gzpo + satisfying the same hypotheses as 

(ii) We use Corollary |4.4| so we need to check the two hypotheses with k - N + 2 and 
h = 1 . As in the corollary, we define for ; e N, 

v!-l i 



n :- 



and t; := t 



Using the conditions on r, zi, and fj, simple computation says that for < rj < 1/2, 

Bl-r N+i {Zl) D B 2R (Zl) D B R (Zo), 



, 16/? 2 2fl 2 16# 2 

f ^ < R + - A — 2 < — < —r 



< -R L + 



2 4(4 + 7] 2 )R 2 



< tN+2- 



Thus we have B 2? ^(z x )x(t N J N+2 ) D B s (z )x[^, ^] D fi s (z )x[y, flforO < r, < \ 
and hence inf u < inf u < 1 . We remark that r^ + 2 is comparable to R. 

Now, it suffices to show for some large M' t] > M n , and small < s' n < s n , we have 



(16) 



mf \rr l 



N+2-i 



VI < i < iV + 2, 



where and are the constants in Corollary |4.4| We notice that Bp x?N+2 (z ) c B ai?N+2 (zi) c 
Si2.i-»(Zo) and 

M - I 2 64/? 2 3 
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since d(z ,Z\) < R,\h\ < R 2 and ^R ^ r N+2 < ^R- Then for i = 1,2, ■• ■ ,N + 2, we have 



-2(n+l) 



|/+| „ + I < 



2(n+l) 



^2(«+l) 



7 11/ 



+ n«+l 



2(n+l) 



(2<)" 



B ?sfe )x(-3/? 2 ,^ + /? 2 ] 



|B4!«fe)l« 2 



< C(n, 77V — < C{n,rf)s' 

n \B ai - n (zi)\n 



l „ + i\ B PCr N+ ,Mo)\r h 



\B ai7i (ziWt 



where we use that f ; R < < ^iR and the volume comparison theorem in the last 

inequality and the constant C(n, rj) > depending only on n and 77, may change from line 
to line. Since d(z„,zi) < R, we use the volume comparison theorem again to obtain 



-2(n+l) 



f ' 



in n+1 <C(n,y7) e ; 



, n+ il B Ar w+2 (z )|r; 



■2 

Af+2 



\B aifi {ziWi 



< C(n,r])E n 



,n+\ \ B "ir N ^SZ\)\ r N+2 



\B ai T.(ZlWf 



(- \ n+2 
r N+2 \ \ r n+1 

-z— I <C(n,rj)s n 



(n+2)(N+2-i) 



We select M' > M n large and < s' < s n small enough to satisfy 



C(n,rj)s' n n 



(«+2)(W+2-i) 



s jti if W + M i VI < i < N + 2, 

which proves ( [16) . Therefore, Corollary 4.4 (after translation in time by t\) gives 

In < M' N+2 \ n K^zutA 

i ' > > 0. 

l^r.^Czi.^i +4r 2 )]| 



5. Parabolic version of the Calderon-Zygmund decomposition 

Throughout this section, we assume that a complete Riemannian manifold M satisfies 
the condition We introduce a parabolic version of the Calderon-Zygmund lemma ( 
Lemma [5T7] ) to prove power decay of super-level sets in Lemma 6.1 (see ifWllCallCClO . 



Christ |Ch| proved that the following theorem holds for so-called "spaces of homogeneous 
type", which is a generalization of Euclidean dyadic decomposition. In harmonic analysis, 
a metric space X is called a space of homogeneous type when X equips a nonnegative Borel 
measure v satisfying the doubling property 

v{B 2R (x)) < AMBr(x)) < +oo, MxeX, R > 0, 

for some constant A; independent of x and R. From Bishop's volume comparison (Lemma 



2.2 1, a complete Riemannian manifold M satisfying the condition ([3]l is a space of homo- 



geneous type with A\ =2™. 
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Theorem 5.1 (Christ). There exist a countable collection \Q k ' a c M : k e Z, a e 4} of 
open subsets of M and positive constants < Sq < 1, C\ and c 2 {with 2c\ < c 2 ) that depend 
only on n, such that 

(i) \M\U a Q kM \ = 0forkeZ, 

(ii) ifl<k,ae h, and B e I,, then either Q k ' a c or Q k - a n Q l # = 0, 
( Hi) for any (k, a) and any I < k, there is a unique B such that Q k,a C Q ^, 

(iv) diam{Q k - a ) < c 2 6 k , 

(v) any Q k,a contains some ball B ClS t(z k ' a ). 

For convenience, we will use the following notation. 

Definition 5.2 (Dyadic cubes on M). 
(i) The open set Q = Q k ' a in Theorem 5.1 is called a dyadic cube of generation k on M. 



From the property (Hi) in Theorem 5.1 for any (k,a), there is a unique fi such that 
Qk, a c Qk-ifi We ca u Qk-hP the predecessor of Q k ' a . When Q := Q k - a , we denote the 
predecessor Q k ~ l '@ by Q for simplicity. 
( ii) For a given R > 0, we define ka G N to satisfy 

c 2 5 k f l <R< c 2 5 k f 2 . 

The number k R means that a dyadic cube of generation k R is comparable to a ball of 
radius R. 

For the rest of the paper, we fix some small numbers; 

^3M(0,«, ft : =*^ e (o,|), 

rj := min((5, 5,) e ^0, M and k := | yjl - 

By using the dyadic decomposition of a manifold M, we have the following decompo- 
sition of M x (7*1, T 2 ] in space and time. For time variable, we take the standard euclidean 
dyadic decomposition. 

Lemma 5.3. There exists a countable collection \K k - a c M x (T\,T 2 ~\ : k € Z,a e J^} 
of subsets of M X (Ti,T 2 ] C M X R and positive constants < <5o < 1. c\ and c 2 (with 
1c\ < c 2 ) that depend only on n, such that 

(i) \M x (Ti , T 2 ]\ \J a K k ' a \ = Ofor keZ, 

(ii) ifl<k,ae J k , and j3 e then either K k - a c K lfi or K k ' a n K lfi = 0, 
( Hi) for any (k, a) and any I < k, there is a unique f3 such that K k ' a C 

(iv) diam(K k ' a ) < c 2 6 k x cjsf, 

(v) any K k ' a contains some cylinder B CiS k(z k ' a ) X (t k ' a - c\8^, t k - a ]. 

Proof. To decompose in time variable, for each k e Z, we select the largest integer Nt e Z 
to satisfy 

4 C 2°0 S 2 2N k < C 2°0 ■ 

For k-th generation, we split the interval (Ti , T 2 ] into 2 2Nk disjoint subintervals which have 
the same length. Then we obtain \ Jk\ = \h\ ■ 2 2Nk disjoint subsets on M x (T\ , T 2 ] satisfying 
properties (i)-(v). □ 

For the rest of this section, let \K k - a C Mx (T 1 ,T 2 ] : k e Z, a € J k ] be the parabolic 



dyadic decomposition of Mx (T t , T 2 ] as in Lemma 5.3 
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Definition 5.4 (Parabolic dyadic cubes ). 

(i) K — K k ' a is called a parabolic dyadic cube of generation k. If K :— K k - a C K k ~ l 'P —: 
K, we say K is the predecessor of K. 

(ii) For a parabolic dyadic cube K of generation k, we define l(k) t o be the length ofK in 
time variable, namely, l(k) — for M X (7\, T 2 ] in Lemma . 



5.3 



We quote the following technical lemma proven by Cabre MCal Lemma 6.5]. 
Lemma 5.5 (Cabre). Let z e M and R > 0. Then we have the following, 
(i) If Q is a dyadic cube of generation k such that 

k> ka and Q c Br(z ), 
then there exist Z\ 6 Q and r^ € (0,R/2) such that 



(17) B Sn (zi) cQcQc B 2n (zi) c Bu ( Zl ) c Bu R ( Zo ) 

I] K '1 

and 

(18) B iR ( Zo ) C Bn R ( Zl ). 

In fact, for k > kn, the above radius rt is defined by 

r k 2 C2 ^o = j s o- 

(ii) If Q is a dyadic cube of generation ka and d(z a , Q) < 8\R, then Q C -Sr(z c ) and hence 
( fT7| ) and ( fT8| l hold for some z\ e Q and r^ e f )• Moreover, 

B SlR (z„) c B 2rtK (zi). 

(Hi) There exists at least one dyadic cube Q of generation kn such that d(z a , Q) < 8\R. 
We remark that for k > kR, 

1 

4 ( 

and ( fTT) gives that for any aeR, 



(19) K m (zi,a) cQx(a- l(k),a] c K 2rk (zi,a) 

Definition 5.6. Let m e N. For any parabolic dyadic cube K :— Q x (a — l(k),a] of 
generation k, the elongation of K along time in m steps (see HKLI ), denoted by K , is 
defined by 

K :- Qx (a, a + m ■ l(k - 1)], 

where l(k) is the length of a parabolic dyadic cube of generation k in time and Q is the 

— m 

predecessor of Q in space. The elongation K is the union of the stacks of parabolic 
dyadic cubes congruent to the predecessor of K. 

Now we have a parabolic version of Calderon-Zygmund lemma. The proof of lemma is 
the same as Euclidean case so we refer to IfWll for the proof. 

Lemma 5.7 (Lemma 3.23, fWl ). Let K\ — Q\x(a- l(ko), a] be a parabolic dyadic cube of 
generation ko in M X (T\, T 2 ], and let < a < 1 and m e N. Let 3\. C K\ be a measurable 
set such that \ J{ C\ K\ \ < a\K\ \ and let 

Jl™ := U {F" : \Kn M\ > a\K\, K, a parabolic dyadic cube in n (Q x x R) . 

Then, we have 

117 

(m+ \)a 
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6. HARNACK INEQUALITY 



In order to prove the parabolic Harnack inequality, we take the approach presented 
in [W| and iterate Lemma |4.3 | with Christ decomposition (Theorem 5.1 1 and Calderon- 



5.7 1. We begin this section with recalling that r\ e (0, |j is 



Zygmund type lemma (Lemma 

fixed as in the previous section. So the uniform constants fi, p s' n and M' n in Proposition 4.6 
are also fixed and we denote them by fi, sq and Mq for simplicity. 
We select an integer m > 1 large enough to satisfy 

m 1 



(m+l)(l-/z) 



1 

1 2 



where < fi < 1 is the constant in Lemma 4.3 For T\ := -3R 2 and T2 '■= (rf + 1 + m^R 2 , 



we consider a parabolic dyadic decomposition of M x (Ti, T2] in Lemma |5T3] and fix the 
decomposition for Section[6] 

6.1. Power decay estimate of super-level sets. 

Lemma 6.1. Suppose that M satisfies the conditions ([3j,ffl. Let z £ M, R > and t e 
[3, 16]. Let u be a nonnegative smooth function such that J£u < f in B5o R (z )x(—3R 2 , . 
Assume f/za? 



B«(z„)x 



inf m < 1 

2ft 2 rfl 2 1 
,2 ' , 2 J 



and 



R l 



B50 s (z o )x(-3fi 2 ,^ 



21 I/T +I 



for a uniform constant < £\ < Eq. £e? K\ be a parabolic dyadic cube of generation k R 
such that 

Ki := Qi x ( fl - l(k R ), t x ] c gj x (-fl 2 ,/? 2 ), 

where Q\ is a dyadic cube of generation k R such that d(z , Q\) < <$iR- Then for i — 1,2, ■ ■ • , 
we have 



(20) 



\{u > M[}r\Ki 



(>-?)' 



where < S\ < eq and M\ > depend only on n, A, A, and a L . 



Proof, (i) As Proposition 4.6 we use Lemma 4.5 to assume that a nonnegative smooth 
function u defined on Bi9 R (z ) x (T\, T{\ satisfies that inf u < 1 and J^fu < f in 

" 3 B«fc,)x[ 

B^ R (z„) x (T\, T2] for some / with 



I 2«2 T«2| 

^ >, 2 



1/ 



+ 1)1+1 



B 49 Jz^x(T,,T 2 ] 



50 

< — £1 < 2ei . 
49 



(ii) According to Lemma 5.5 there exists a dyadic cube <2i c B R (z ) of generation 
such that d(z ,gi) < 6^R. We find zi € Q x and r fe e [ffl, ±i?) satisfying (|T7]> , (fl~8]> 
and B SlR (z„) c ^(zi)- Since 77^ < l(k R ) < 4r 2 kR+l = 45 2 r 2 s < (5 2 /? 2 , we find h e 
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(-R 2 + l(k R ),R 2 ) such that K\ := Q\ x(t]-l(k R ), t\ \ is a parabolic dyadic cube of generation 
k R of M x (Ti , T 2 ]. From ([19]), we also have that 



K nrkR {Zl,h) C^C K 2rkR {Zl,h). 

We use the induction to prove pO) so we first check the case i = 1. We notice that 

d( Zo ,z 1 )<R,r kR e [fK,itf)c and |fi| < fl 2 . Weset £l := ( ^ff^ )^ f . Then, 

m satisfies the hypotheses of Proposition 4.6 with r = and AT = 1, so we deduce that 

\{u<M 3 }nK m (zuh)\ \{u<M 3 }nK nr (zi,h)\ \{u < Mh n kA 

< f 1 ^ 77, ; : TTT = TP 7 7\i < 



l^,«y»Czi,*i+4^)| \K airkR , air 2(zuh)\ ' |*i 



Thus, we have for M\ > M; 



k R " 1 air kR ,a 2 q R ^ 

0' 

|{K>Mi}n#i| 

< 1 — u < 1 — . 

I*i I 2 

(iii) Now, suppose that p0] ) is true for /, that is, 

|{w > M\) r\K 



H)' 



l*i I 

To show the (i+l)-th step, define for h > 0, 

S h ■.= {u>h}nB« LR (.Zo)x(TuT 2 \. 



Pm; n * / u\i 

We know < 1 . If h > is a constant such 

l*i I V 2^ 



that 



l*i 

^>(l-^)' +1 for yi:=S /lM , n/v ; . 



then we will show that h < M\ for a uniform constant M\ > Mq > 1, that will be fixed 
later. 

Suppose on the contrary that h > M\. From (ii), we have ■ — - < ■ — ■ < 1 - fx for 

l*i I l*i I 

Mi > M ( , and h > 1. Applying Lemma 5.7 to J[ with a — 1 - ju, it follows that 

1 



11 (m+l)(l-ju) 1-f 



We claim that 



(21) ^ c 

for /z > C\M™ > 1, where a uniform constant C\ > will be chosen. If not, there is a point 
(jcj, s\) e ' \£ui> an< i we find a parabolic dyadic cube K := Q x (a - l(k),a] c of 

generation £(> A^) such that 

\ftC\K\>(l-fi)\K\ and (jq, si) e 



from the definition of ^t'"^- According to Lemma 5.5 there exist z\ e 2 c 6i c B R (z ) 
and e (0,R/2) satisfying ((TTJ, ([18]), K m {z\,a) c K c K lrk (z\,a) and 

e A = Qx(a,a + m ■ l(k - 1)] c B2r t (Zi) x(a,a + m- Ar\\. 



24 



SEICK KIM, SOOJUNG KIM, AND KI-AHM LEE 



We note that 



inf u < u(x\, s\) < 



B 2 r k (ziMa,a+m-4r 2 k ] 



hM\ 



and 



B airk (zi) X («- fa 2 + ifl*)r 2 k ,a + m ■ 4r t 2 ] c B aiR { Zo ) x (-3^,(1 + m)R 2 ] , 
since r k < R/2 and a e (ti - l(k R ), t\] c (-R 2 ,R 2 ). We also have that for j = 1, • • • , m, 



(22) 



2 (z, 



,a+(m— /fl)-4f?) 



ir 



< £0" 



m: 



m—j+1 




Indeed, the volume comparison theorem and the property ( fT8) i will give that 



\l 



\f 



+ i«+l 



•k 'k 



RSI .R'n 



R 2 

\Bp lR (z )\^ (a 2 R 2 y» 
CiR 2 /2 



"■^ "L" +1 (Bi2 K (z„)x(r,,r,]) 

urn 



-iiniL-i 



B49(&)x(r 1 ,r 2 ] 



fii?fc)x(riJ 2 ] 



-urn 



S4 B49 fe,)x(ri,r 2 ] 



< Ciei, 



where a uniform constant Ci > 1 depends only on 77,71 and m. For h > C\M™ and M\ > 1, 
we have that 



1 ^ 2 



ir 



2 (zi,a+(m-j+l)-4r z k ) 



hM\ hM\ 
< Ciei < s -—- < e 



M m - m; 



"1-7+ 1 ' 




which proves ( |22] >. Thus, we can apply Lemma 4.3 iteratively to My := ° Mi u, for 1 < j < 
m, to deduce 

\\u < /zMj} n ^ rt (zi,a)| \[u < /iMj} n A"! 



\K airk , air l(zi,a + 4r 2 k )\ 



\K\ 



However, this contradicts to the fact that ffl n K\ > (1 - fi)\K\. Therefore, we have proved 



that flf_ M c £ hM , for h>C\M™. 
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(iv) Since IS^O^I < (l - f) ! |^i|,wehavethat|S w ; n *i| < |fl*jn*i| < (l - |J |*i| < 



r\9{\ for ft > CiM™. Then, by using ((21}, we obtain 



to 



(m+l)(l-A0 
m 



1^1- 



S„„ ; OA", 



with a- := 



m _ 1 
(m+l)(l-/j) 1-g 



> 0. We find a point {x\, s\) e 3[™ ^\K\ and a parabolic dyadic 



cube := Q x (a - Z(fc),a] c of generation k(> k R ) such that (xi, si) e K , and 
n K\ > (1 - n)\K\. We may assume that 

*1 > *1 + " ' 

since c gi x ( f i _ + 00 ) an d 



— 'jq^I ' > a(l - j J Using Lemma 



5.5 



again, there exist Z] e Q c c B R (z ) and e (0,R/2) satisfying (fT7]i,([T8j), ana 

K nrk (z\,a) c K c K2 rk (zi,a). Then we have 



V8^V 2/ 4a/2^\ 2/ 12/ 



si < a + m • l(k -l)<?i+ to • 4rt 



and hence 



for a uniform integer N > independent of i e N. We apply Proposition 4.6 to u in order 
to get 



|{u < Mf+ 2 } n K Vk (zi , a)\ \{u < M ( Q N+2)i } n K\ 



\K\ 



sincere > (if' R, and(zi,a) e K c K { c B R (z M-R 2 ,R 2 ). Ifh > Mi := max{CiM2',M^ +2 }, 
this implies 



1 -p. > 



[u > M { N+2)i ] n id |{m > fcMj } n k\ \ji nK \ 



\K\ - \K\ \K\ 

which is a contradiction to the fact that \M n K\ > (1 - fi)\K\. Thus, we have h < M\ for 



a uniform constant M\ := max{CiMg", M^ +2 }. Therefore, we conclude that 
(l — f ) , completing the proof. 



\ju>M^)nKi\ 

ft I 



The following corollary is a direct consequence of Lemma 6.1 which estimates the 
distribution function of u. 

Corollary 6.2. Under the same assumption as Lemma |oT7] we have 

\{u > h}nKi\ 



(23) 



1*1 



< dh~ e \/h > 0, 



where d > and < e < 1 depend only on n, A, A, ant/ a £ . 



26 



SEICK KIM, SOOJUNG KIM, AND KI-AHM LEE 



Another consequence of Lemma 6.1 is a weak Harnack inequality for nonnegative su- 
persolutions to u = /. 

Corollary 6.3. Under the same assumption as Lemma 6.1 we have for p„ :— |, 
(24) | f w p °| <C, 

\\K KR (ZoM Jk M I 
where C > depends only on n, A, A, and ai. 

Proof. Let k = k R and let := Q k > a x (**•<* - l(k),t k ' a ])^ r be a family of parabolic 

dyadic cubes intersecting K kR (z , 0). For a e J' k , we have that K k ' a c B( K + Sa ) R (z )y-(-R 2 , R 2 ] 
since d(z , < kR(< 6iK), diam(Q k - a ) < c 2 6 k < 6 R, and -R 2 + l(k) < -k 2 R 2 < t k ' a < 
l(k) < 6 2 R 2 . Since 

* (^J Z ■ ^ >- {^hJ Z ■ 

x are/-; 

2 ib^^c^)! • * 2 fl 2 



k + #o 2(^o + 2/c) j 

^•^r^)"Z^)i-^ 2 ' 

v /e + £ 2((5 + 2k)/ ^ 

the number of parabolic dyadic cubes intersecting K kR (z o ,0) is uniformly bounded. 
Thus for some K k - a with a e J',, we have 

f u p ° < \J' k \ ■ f u p ° 

< \J' k \ • + Po J hP°~\u > h] n x^idftj 

< 141 ■ + Po d\K k n J . 

from Corollary |6.2| where d and e are the constants in Corollary |6.2| 
By using the volume comparison theorem, we conclude that 

1 f Po<c \K k ' a \ < (k + 6q\" &i 
\K KR {zoM J^fo,o) " °l^fc,0)|- °\ k ) k 2 

for Co := • jl + p d r 1 "^ 2 ^} since X*» c B ilc+lh)R (z„) x (f^ ff - <5 2 /? 2 , «**]. □ 
6.2. Proof of Harnack Inequality. So far, we have dealt with nonnegative supersolutions. 



Now, we consider a nonnegative solution u of Jz? m = /. We apply Corollary 6.2 as in MCal 
(see also IfWl ) to solutions C\ - C 2 u for some constants C\ and C 2 . 

Lemma 6.4. Suppose that M satisfies the conditions (|3J,Q. Lef z e M, > anof t e 
[3, 16]. Lef u be a nonnegative smooth function such that S£u = f in Z?5o R (z )x(— 3/? 2 , . 
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Assume that 



Br(zM 



inf u < 1 and 

2R 2 jR 2 1 



R z 



fl50 R (z„)xl-3R 2 ,^| 
V if 



* 4 =:8 



/or a uniform constant < e\ < \ as in Lemma ]6J\ 

Then there exist constants o~ > one/ Mq > 1 depending on n, A, A and ai such that for 



v :— ^ M °^ 2 > 1, the following holds: 

If j > 1 is an integer and z\ € M and t\ € R satisfy 

d(z ,Zi)<i<R, \h\<K 2 R 2 

and 

u{z\,h) > v ! ~ x M , 

then 

(0 K fLl ^)^^ B ^o)><{-3R\fl 
(ii) sup u > v ] Mq, 

*>f?H (z " ri) 



where Lj :— crM {) " +2 v "* 2 R andO < e < 1 as in Corollary 
Proof. We select cr > and Mq > 1 large so that 



6.2 



and 



cr> ^-(2d2 e )^ 2 

C\0() 



o-M^- +dMj<-, 



where d, e, c\,ci and 6q are the constants in Corollary 6.2 and Theorem 5.1 Since Lj < 



kR , T]R 



< f , d(z ,zi) <kR<R and hi < k z R z < 



2»2 . T«- 



we have 



B §Lj ( Zl )x\ h -\3 + - 2 )L^ 



tR 

cBso R (z )x\-3R 2 ,^- 

■ TJ 



so (i) is true. 

Now, suppose on the contrary that 



sup u < v j Mq. 



>f?H (z,A) 



Let kj : = kz,. > kR with Lj in Definition 



5.2 



From Lemma 



5.5 



there exists a dyadic cube 



<2z,j of generation kj such that d(z\,Qt) ^ S\Lj. We also find a parabolic dyadic cube K L . 
of generation kj such that 



tL 2 tL 2 



since l(kj) < d^Lj. Let K\ be the unique predecessor of Ki. of generation k^, that is, 

K L . cKr-=Qi x{a-l(k R ),a]. 
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Then we have 

d{z , Qi) < d(z , Q Lj ) < d(z ,zi) + d(zi,Q Lj ) <kR + SxLj < 6 { R 
and (a-l(k R ),a] c (-R 2 ,R 2 ) 

since 

tL 2 16L 2 
l(k R ) + \ti\ + ^-+L 2 < l(k R ) + \h\ + ^+L 2 



Now, we apply Corollary |6.2| to u with K\ to obtain 



16 



16 



64 



R 2 < R 2 . 



(25) 



u > v 



jMo 



riK, 



u > v- 



<d\v 



■Mi 



\K 



il- 



On the other hand, we consider the function 

(v-l)Mo 



w := 



which is nonnegative and satisfies 



/ 



from the assumption. We also have w{z\,t\) < 1 and 

l/l „ l/l 2(M„- 1/2)|/| 



V j- l (v-l)M (v-l)M 



2|/|- 



By using the volume comparison theorem with <^R< \ and Bu 4 ?«(z ) c B50 i«(zi) 

4 8 >? ,2 S ,3 S 

we get 



12/r 

2(7?fl/8) 2 



2L* 



|B» L ,(z.)l-{(3 + ^)L^ 
-ll/lli- 



-II/IIl- 



|fi».^(zi)l^{(3 + ^)(^/8) 2 p 

- — ^ — xll/ll- 

Ifiii4^fe)l^{(3+^)(77W}" +1 
|fin 4s (z )|-{(3 + ^)^}-' l/llL "K B ^ ( < 3S2 '7 



< 2 



1 17- 

50 \^ <V; 

44 



<ei- 



Applying Corollary 6.2 to w in /f L; , we deduce that \{w > M ) n | < dM e |/r Lj |, i.e., 



m < v 
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Putting together with ( |25| ), we obtain 

\K Lj \<2d2^M E \K l \ 

since dM n e < f < 1/2. From TheorempTI there is a point z* e 2l sucn that B i,(z») c 
_ 2 ' cia 



Q Lj <z Q x <zB c ^ s (z,). Then we have 



ciS' 



■ Kkj) < \K L .\ 

< 2d2 e v- Je M Q e \Ki \ = 2d2 e v- je M e \Qi \ ■ l(k R ) 
<2d2^M E \B C2 ^(z,)\-c 2 2 6l kR 



< 2d2 e v- ie M Q t 



from the volume comparison theorem. This means 



C2<V 



B M 



r 2 s 2k R 
c 2°0 



6 k Q ' < (2d2 e )^- M Q " ll v-£i —6 k Q R . 

^* 1 



Since c 2 6 k f 1 < R < c 2 8 K ( « \ we deduce that 



fe-2 



/.; < <'mV* : < -^L (2,/2') • M p - r ^ 



c v 5 l 



o 



C'2 



< (2d2 e ) ^ M Q 11+2 v~—2R< o-M Q " +2 v~^-R = Lj, 
in contradiction to the definition of Lj, Therefore, (ii) is true. □ 
Thus we deduce the following lemma from Lemma [674] 

Lemma 6.5. Suppose that M satisfies the conditions (13), (Hi. Let z 6 M,R > and r € 
[3, 16]. Let ube a nonnegative smooth function such that J£u = f in B50 /? (z o )x( — 3R 2 , . 
Assume that 

inf m < 1 



and 



R- 



for a uniform constant < s < 1 in Lemma]6jt\ Then 



< E 



where C > depends only on n, A, A and cLl, 
Proof. We take j„ e N such that 



sup u < C, 



f,50 r kR , fW 16\ 2 
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We claim that 



sup u < v-i" Mo with Mq > 1 as in Lemma 



6.4 



If it does 



not hold, then there is a point (Zj ,tj ) e B«r(z ) x ^j") such that u{Zj a ,tj o ) > 

yj°~ l Mo. Applying Lemma 6.4 with (zi,h) = (z^, f /o ), we can find a point (zj +i,tj +\) e 



^» L , ; ( 3+ ^) L? ^'^) such that 

u{Zj 0+ i,tj 0+ i) > v J °M . 
According to the choice of j a , we have 



kR kR 

d(z ,Zj o+ \) < d(z ,ZjJ + d(zj„,Zj„ + i) < — + — = kR 



2 2 



and 



K R K R~ 



2 D 2 



l^+il ^ If/. I + h - < — + — < k R 



Thus we iterate this argument to obtain a sequence of points (zj, £,-) for j > j„ satisfying 



d(z ,Zj) < kR, \tj\ < k z R z and u(zj,tj) > v J 1 M , 

CO DO _ „ CO 

ZkR 50 \ 

d(Zi,Z M ) < — + 2j —U < kR and |f,| < \t jo \+ ^ \t~t i+ i\ < 

1=76 l=Jc t=J 

K 2 R 2 V"1 / T \ 2 2 2 

— h ^ |3 H — j jLj < k'R for j > j . This contradicts to the continuity of u and 



therefore we conclude that 



sup u < v ! ° 1 Mq. 



Now the Harnack inequality in Theorem 6.6 follows easily from Lemma 6.5 by using a 
standard covering argument and the volume comparison theorem. 

Theorem 6.6 (Harnack Inequality). Suppose that M satisfies the conditions dJhfi}. Let 
Z £ M, and R > 0. Let u be a nonnegative smooth function in K2r(Q,4-R 2 ) C M X R. Then 



(26) 



sup u<C\ inf u+R 2 \\ \&u\ n * l \ \ , 

K R (z ,2R 2 ) U*(z„,4« 2 ) \Jk 2 r(z,„4R 2 ) I I 



where C > is a constant depending only on n, A, A and a L . 



Proof. According to Lemma 6.5 for t e [3, 16], a nonnegative smooth function v in 
K,„ i \ , (x, f + ^ ) satisfies 



(27) 



since | < 1. 



sup v < C 



inf v + r 



Now, let (x,t) e K R (z„,2R 2 ) = fi«(z ) x (fi 2 ,2fl 2 ] and (y,s) e tf«(z ,4fl 2 ) = fl R (z ) x 
(3R 2 ,4R 2 ]- We show that 

R 2 



u(x, t) < C \ u(y, s) + 



\K 2R (z ,4R 2 )\ 



— \\J£ «llL" +1 (/f2«(z ,4R 2 )) 
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for a uniform constant C > depending only on n, A, A and «£. We consider a piecewise 
C 1 path y : [0,1] — > M, y(Q) = x,y{l) — y, I < 2R, consisting of a minimal geodesic 
parametrized by arc length joining x and z , followed by a minimal geodesic parametrized 
by arc length joining z and y. We notice that y([0, I]) c B R (z c ,) and d(j(si), y(s 2 )) < 
\si - s 2 \. 

We can select uniform constants A > and iVeN such that 



A := max 



16 

, AT, we define 



(*,-,*,-) := \y\i 



s - t 
~N~ 



+ t\eB R (z„)x[R 2 AR 2 ]. 



Then we have (xq, to) = (x, t), (x N , t N ) = (y, s) and for i = 0, ■ ■ • ,N — 1, 

s I 2R (A kR kR 
d{Xi+\,Xi) < — < — < — z — ■ — < — , 
V N N 3kt] 2 A 2A 2 A 

3R 2 R 2 _s-t 3R 2 16R 2 

-^<-< t M -?,- — < — < 

We also have that K 50R # (xtJi) c K 2 r(z AR 2 ) for i = 1, • • • ,N since < /?. 

We apply the estimate ( |27] > with r = f,T = (f,-+i - fj)^r- and (jc, f) = (Xi+i.fj+i) for 
i = 0, 1, • • • , N — 1 and use the volume comparison theorem to have 



u(xi, ti) < C 



< C 



u(x M ,t M ) + 



(R/A) 2 



u(x i+ \,t i+ \) + 

< C<u(XM,t M ) + 



\K 50B ff 2 (xi + i , ti+{}\ » 

(R/A) 2 
B rx (x M ) -(3 + ^)1^ 

R 2 



-\\^u\\l"*'(K 2 r(ZoAR 2 )) 



-||-5f"llL» +1 (^2Kbo,4J? 2 )) 



|B 3 *(x, +1 )-4* 2 |; 



ll-^ M llL» +1 (A: 2s fc,4fi 2 )) f . 



where a uniform constant C > may change from line to line. Since Bj^Xi+i) D B2r(z ), 
we deduce that 



Therefore, we conclude that 

u(x, t ) < C < u(y, s) + 



R 2 



\B 2 r(Zo) ■ 4R 2 \" 



-\\^u\\l"+ 1 (k 2 r(z AR 2 )) 



R 2 



\K 2R (z ,4R 2 )\^ 
for a uniform constant C > since N eN is uniform. 



■\\-^u\\l"+KK2r(z AR 2 )) 



Arguing in a similar way as Theorem 6.6 Corollary 6.3 gives the following weak Har- 
nack inequality. 

Theorem 6.7 (Weak Hamack Inequality). Suppose that M satisfies the conditions ([3]),(|4|. 

Let z 6 M, and R > 0. Let u be a nonnegative smooth function such that I£u < f in 
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K 2 r{ZoAR 2 \ Then 



f u p°Y <c{ inf «+Wf \rr'Y 

JK R (z„,2R 2 ) I \K R (z,,AR 2 ) \Jk 2r (z„AR 2 ) I 



where < p a < 1 and C > depend only on n, A, A and aj_.. 



Proof. Let e > be the constant in Corollary |fT2| and let p (1 :- |. We consider a parabolic 

Let k :- kuR for the constant 



decomposition of M x (0, 4R 2 ] according to Lemma 
A > in the proof of Theorem 



6.6 



5.3 



Let [K k ' a := Q k - a x (t k '° - l(k), t k ' a ]} aeJ , be a family of 

parabolic dyadic cubes intersecting K R (z , 2R 2 ). We note that diam(Q k ' a ) < C2<5q < So ■ ^ 
and < £g • '^rr- Following the same argument as Corollary 6.3 we deduce that \J' k \ is 
uniformly bounded and 



I 



u p ° < \J\ 



K R (z,„2R 2 ) 



for some K k - a with a e J' k . Then we find (x, t) e K k > a n B«(z c ) x [/J 2 , (2 + 5 2 £)R 2 ] such 
that c K&ix, t) since diam(Q k ' a ) < S Q ■ f and Z(/fc) < S 2 ■ Since <i(z„, x) < and 
B*r(x) c fi/' 1+ £\o(z ), we have 



(28) 



— 5- f u* < — ^— [ u"° 

\K R (Zo, 2R 2 )\ Jk^^) \Km (x, t)\ J KkR m 

A x 



for Co := \J' k \ (l + j) ■ by using the volume comparison theorem. 
We set 

inf it —: u(y, s) 

K R (z„AR 2 ) 



for some (y, s) e K R (z AR 2 )- As in the proof of Theorem 6.6 we take a piecewise geodesic 
path y connecting x to y. Let N e Nbe the constant in Theorem 6.6 For i = 0, 1,- • ■ , ./V, 
we define 



{Xi,td:=\y\i-\,i 



s - t 

~N~ 



+ t\eB R (z„)x[R l AR 1 ]. 



Then we have (xq, to) = (x, t), (x^, t^) = (y, s) and for i = 0, 



d{x M , Xi )<--- and 



< U + \ - tt < 



,N-1, 
16R 2 



T] 2 A 2 



It is easy to check that for any i = 0, 1, • ■ • ,N - 1, B^(Xi) n Z?^(x, + i) D fi«s(x,+i) and 
hence 



inf u < inf m < 



(29) 



r 

2^1 l - r 
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On the other hand, Corollary |6.3| says that for i = 0, 1, • ■ • , N — 1 



\Kdt ( 



^ f A 

Xi,ti)\ jKateA) 



l/Po 



< C 



inf u + 

KkR taA+i) 



(fl/A) 2 



^50 « iff 2 , (-X^fi+l) 



1 ||/ + |L"+1(^ S (Z„ 



4R 2 )) 



< C { inf m + 

K»s (x„t M ) 



R 1 



\k 2 r(z AR 2 ) 



II/1L 



L"*HK 2 r(z„AR 2 )) I 



by using the same argument as Theorem 6.6 with the volume comparison theorem. Com- 
bining with ( |29] l, we deduce 

( \ VP* 



\K 



— f 



K«* (a-,0 



R 2 



<C< inf « h — 

(* f ( Ww ) \K 2R ( Zo AR 2 )\^' 



<C\ inf u + 

Km (x N ,t N ) 



< C < u(y, s) + 



R 2 



K 2R (zoAR 2 )\ l 
R 2 



\k 2 r(z AR 2 )\^ 



' \\f IL» +1 (fc ff (z ,4« 2 ))j 

in ! 

\L'<+HK2R{z a AR 2 )) | ' 



for a uniform constant C > since N e N is uniform. Therefore, the result follows from 
d28b. □ 
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